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Abstract. - Tunneling spectroscopy of epitaxial (100) oriented YBa2Cu307_<s films was 
performed using an STM at 4.2 K. On atomically smooth areas, tunneling spectra revealing 
clear U-shaped gaps with relatively low zero bias conductance were measured. These spectra 
can be well fitted to the tunneling theory into a d-wave superconductor only when introducing 
a strong dependence of the tunneling probability on the wave- vector k. Possible origins for this 
fc-selectivity in STM measurements will be discussed. On other areas, V-shaped gaps as well as 
zero bias conductance peaks are observed, indicating relaxation of fc-selectivity and the effect 
of nanofaceting, respectively. 



The d-wave symmetry of the order parameter of YBa2Cri307_5 (YBCO) has a 

special signature in the tunneling spectra measured for these superconductors. The hallmark 
of this is the zero bias conductance peak (ZBCP) that appears when tunneling along or near 
the nodal direction ([110] for YBCO). This reflects the existence of Andreev bound states at the 
Fermi level residing on the pair-breaking (110) surface |4I5I6I7I<5] . The ZBCP was measured by 
numerous groups, using both macroscopic and microscopic (STM) tunnel junctions 9, lUl lllI 
H2U13lH^lll5| . Manifestation of the d-wave order parameter in the tunneling spectra acquired 
along the anti-nodal ([100] for YBCO) and the c-axis directions is more controversial. The 
theory of tunneling into a d-wave superconductor predicts a V-shaped gap structure for both 
tunneling directions in the case where the density of states (DOS) is isotropically averaged 
over the Fermi surface JSJE]. Here, quasi-particle excitations are observed at any finite energy, 
since the nodal directions in fc-space are monitored, resulting in a linear dependence of the 
dl/dV vs. V tunneling spectra on energy around the Fermi energy (zero bias). Such isotropic 
tunneling scenario is believed to apply to STM measurements |16| . In contrast, U-shaped 
gaps are predicted to appear when fc-vectors along the anti-nodal direction are preferentially 
monitored in the tunneling process |7lllllll4lll7lll8| . Resolving this issue is thus important for 
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understanding STM spectroscopy in high temperature superconductors and other anisotropic 
systems. 

Experimentally, V-shaped gaps are indeed commonly observed in c-axis tunneling measure- 
ments [ZlinilllliniEniEIlES- Recently, however, U-shaped gaps were observed in tunneling 
spectra acquired on the (001) Bi2Sr2CaCu20s+5 (BSCCO) surface ^Jj. This was accounted 
for by a 'fc-selection mechanism' inducing preferential tunneling along the anti-nodal direc- 
tions, resulting from the large overlap, and consequently strong coupling, between the tip's 
electronic states and the 3d orbitals of the Cu atoms, through the mediation of their 4s 
states [23]. However, a corresponding effect was not yet established for the o-axis surface, 
where direct coupling to the Cu atoms is possible. 

The data reported for a-axis tunneling is much more diverse, possibly due to problems 
in preparing well-defined (100) cuprate surfaces, and keeping the exposed Cu-0 planes from 
deterioration. A wide range of features were observed, including zero bias peak structures, 
mostly inside gaps {11113124] . as well as V-shaped |25I26| and U-shaped gaps |18I27| . The zero 
bias conductance peak may result from nano-faceting or roughness of the measured surface, as 
has been demonstrated by various groups .14,21,28,29 . The U-shaped gaps were accounted 
for by assuming a narrow tunneling cone centered along the anti-nodal direction, although, 
as will be shown below, the 'orbital coupling' model discussed above may apply as well. 

In this paper we present an extensive scanning tunneling spectroscopy study of (100) 
YBCO surfaces, directed at resolving some of the questions arising from the contradictory 
experimental findings described above. We have measured films prepared using two different 
deposition techniques, DC sputtering and laser ablation. For both film types, spectra exhibit- 
ing U-shaped gaps with relatively low zero bias conductance are found on smooth areas. These 
spectra are well fitted to the tunneling theory into a d-wave superconductor only when in- 
troducing a tunneling probability that is preferentially strong along the anti-nodal directions. 
However, we cannot distinguish between the 'orbital coupling' and the 'tunneling cone' mod- 
els. In areas showing a rougher topography, only V-shaped gaps and zero bias conductance 
peaks (ZBCP) of varying strength are observed. 

Epitaxial thin YBCO films of various types were prepared by either DC sputtering or laser 
ablation, with relatively broad transitions, onset at 86 (88) K and ending at 82 (81) K for 
the DC-sputtered (laser ablated) films. Both preparation methods produced films yielding 
very similar spectroscopic data in spite of the different morphology. YBa2Cu307_6 with 5% 
Ca doping were grown on LaA103 substrates by DC off-axis sputtering, applying the method 
described in Ref. 30 for c-axis films. Here, however, the deposition temperature was reduced 
to 720° C in order to promote the development of a-axis crystallite outgrowths [31]. These 
films consisted of dense arrays of rectangular a-axis crystallites, typically 20 nm in height 
and a few hundreds of nm in size, covering 70-90% of the surface area, as confirmed by X-ray 
diffraction, Atomic Force Microscopy (AFM) and Scanning Electron Microscopy (SEM). Many 
of these crystallites had large atomically smooth areas on the top. The STM and AFM images 
in Figs, ^ and^ 3 exhibit the surface morphology of these films, showing a single crystal (and 
a junction with its neighbor) and an array structure, respectively. We have also measured 
DC sputtered a-axis YBCO films prepared using the procedure described in Ref. |23|, also 
showing large flat areas. 

The laser-ablated films were grown on (100) SrTiOs wafers by two different dual-step 
processes. First, a 30 nm thick YBCO layer was prepared at a wafer temperature of 600°C 
(or 650°C). Then, a second 75 nm thick YBCO layer was deposited at 750°C (or 785°C). 
The films prepared at the higher temperatures displayed a-axis crystallite outgrowths on a 
c-axis background, similar to the DC sputtered films, as depicted in the STM image of Fig. 
[Ifr. The lower deposition temperature produced films with two coexisting a-axis phases on 
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Fig. 1 - STM and AFM images portraying the surface morphology of YBCO films containing various 
degrees of c-axis (background) and a-axis (outgrowth) orientations, (a) and (b) - DC sputtered YBCO 
films with about 70% a-axis orientation, where (a) is a 180 x 180 nm STM topographic image focusing 
on a single a-axis crystallite, 25 nm in height, and (b) is a 1.6 x 1.6 /im 2 AFM image showing an array 
of outgrowing crystallites, having atomically smooth areas on the scale of a few tens of nm. (c) and 
(d) - laser ablated films, where (c) is a 35 x 35 nm STM image displaying a single outgrowing a-axis 
crystallite, 20 nm in height, and (d) is a 2 x 2 fim 2 AFM image of a 95% a-axis film exhibiting two 
types of areas. The first is composed of small crystallites a few unit cells in height, while the second 
consists of large areas, atomically smooth on scales of 100 nm. 

about 95% of the film area (as determined by X-ray diffraction). One phase consisted of small 
crystallites, a few unit-cells in height, while the other is composed of large areas, atomically 
smooth on scales of hundreds of nanometers, as depicted by the AFM image in Fig. ^i. 

The samples were transferred directly from the growth chamber into a chamber that was 
over-pressured by dry oxygen, then mounted within a few hours into a home made cryogenic 
STM and cooled down to 4.2 K in He exchange gas, being exposed to air for less than 15 
minutes. A total of 5 sputtered films and 4 laser-ablated films were measured. The tunneling 
dl/dV vs. V curves were acquired either directly by the use of conventional lock-in technique, 
or by numerical differentiation of the measured I-V curves, with similar results obtained 
by both methods. We have confirmed that the measured gaps and ZBCP features were 
independent of the STM voltage and current setting (before disconnecting momentarily the 
feedback circuit). This rules out the possibility that the gap features are due to the Coulomb 
blockade, which is known to be sensitive to these parameters All the STM data presented 
in this paper were acquired at 4.2 K with (normal) tunneling resistances between 100 Mf2 to 
lGfi. 

On all the measured samples we have found atomically smooth areas, on the scale of 
tens of nanometers or more, where the dl/dV vs. V curves featured mainly U-shaped gaps, 
with relatively low zero bias conductance. The gap value distribution in these regions was 
relatively narrow, varying between 16 to 18 meV, for both laser ablated and DC sputtered 
films. Typical measurements are presented in fig. [2] (empty circles), where spectrum |2K was 
measured on the laser-ablated film, on the a-axis crystallite shown in fig. \I}p, and spectrum 
[2t> is a representative example for the sputtered films, acquired on the crystallite shown in 
image QJi. The inset of fig. Et presents the I(V) curve from which the differential conductance 
was derived, depicting a BCS-like structure. The red curves were calculated using the theory 
for tunneling into d-wave superconductors |HJ|7], taking an equal weight for all fe-vectors, 
i.e., an isotropic tunneling matrix. It is clear that the V-shaped gap structure thus obtained 
deviates considerably from the experimental curves at the low energy regime, the experimental 
curves being by far more flat. The blue curves in fig. [21 on the other hand, were calculated 
assuming preferential tunneling in the anti-nodal directions. In fig. |2K we used the 'orbital 




Fig. 2 - a) A tunneling spectrum measured on the a-axis crystallite shown fig. ^ ( em Pty circles), 
with the corresponding I-V curve shown in the inset. The blue curve is a fit using the 'orbital 
coupling' model with a 16meV gap and V — 0.06A. b) A spectrum measured on the crystallite shown 
in fig. [Ql (empty circles), fitted to the 'narrow tunneling cone ' model, A = 17meV and T — 0.1 A. 
An equally good fit was obtained for the 'orbital coupling' model, shown in the inset. For comparison, 
we present spectra calculated assuming isotropic averaging over the Fermi surface, with the same A 
and r values, showing a large deviation from the measured spectra at the low bias region (red curves). 



coupling' model, whereas in fig. |5Jd the 'narrow tunneling cone' approach was applied, both 
reproducing well the measured U-shaped gap structure. 

For the 'orbital coupling' model we used the tunneling matrix element given by Misra et 
al. J3> \M(k)\ 2 oc \cos(k x a) — cos(k y a)), where k x and k y are the wave vector components 
along the main axes in the a-b plane. Note that this model was applied in Ref. [T7J to c-axis 
tunneling. However, it should also be valid for tunneling into the (100) surface, where the 
lobe of the Cu ci-orbital is perpendicular to the surface, thus strongly overlapping directly 
with the tip states. The blue curve in fig. [2^ was calculated using this model, with a gap 
of A = 16meV and a small Dynes broadening parameter V — 0.06A, showing good 
agreement with the experimental data in the gap region. Equally good fits were obtained 
using the 'narrow tunneling cone', described below. 

The 'narrow tunneling cone' model is, in a sense, a semi-classical approach. Within the 
WKB approximation, and at bias voltages much smaller than the tunneling barrier height $, 
the tunneling transition probability decays exponentially with increasing fey, the transverse 
component of the wave vector k: p{k) oc exp[—{h 2 d/(2m^) 1 ^ 2 }k 2 -,], and consequently p{9) cx 
exp[— (3 sin 2 (9)]. Here, 8 is the angle between the momentum of the tunneling electron and 

the surface normal and /3 = \J 2m/h 2 , where Ep is the Fermi energy and d is the width 
of the tunneling barrier ^2E >34 35 . For typical tunneling junctions, this will result in a 
narrow tunneling- momentum cone around the surface normal (A6> ~ 20°, full angle). 

This type of reasoning doesn't alter the characteristics of tunneling in the [110] or the [001] 
direction, but predicts a U-shaped gap for the (100) surface. Figure 03 demonstrates a fit 
obtained using this model with A = 17meV and T = 0.1A. This particular fit was calculated 
with A9 — 20°, but no significant difference is found in the range of 0° < A8 < 30°. In 
the inset we present, for comparison, a nearly identical spectrum calculated using the 'orbital 
coupling' model. We were not able, in this work, to differentiate between these two models, 
but it is clear that a fc-selection mechanism is essential to account for our results. 

It should be pointed out that V-shaped curves obtained from the isotropic tunneling model, 
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Fig. 3 - Tunneling spectra displaying V-shaped gaps, a) Spectrum measured on a rough region of 
the laser ablated sample shown in fig. [fll. Inset: Spectrum measured on an a-axis outgrowth, b) 
Spectrum taken on a sputtered film. In both a) and b), the dotted curves are fits calculated assuming 
an isotropic tunneling probability, with gap values indicated in the figure and V = 0.1 A for the inset 
and about 0.2A for the other two curves. 



such as those presented in red in fig. [21 can be smeared using very large Dynes broadening 
parameter until a structure resembling U-shaped gaps would emerge. Then, however, the zero 
bias conductance becomes much larger than that of the experimental data. 

In areas having a larger degree of surface roughness and disorder, U-shaped gaps were never 
observed, and the spectra showed only V-shaped gaps or ZBCPs. The zero bias conductance 
in these curves were typically larger as compared to the U-shaped spectra, varying between 
20% to 50% of the normal differential conductance (at high bias). Best fits of these data are 
obtained using the isotropic tunncling-probability model with a lifetime broadening of around 
0.2A. Fits with any of the fc-selective models described above, yield unreasonably large Dynes 
parameters of more than 0.5A. In fig. |2|we present typical V-shaped gap structures measured 
on various film types, along with the corresponding fits depicted by the dotted curves, with 
gap values denoted in the figure. The spectra in fig. were measured on the two types of 
laser-ablated films. The one in the mainframe was taken on a rough part of the film shown 
in fig. nji, and the other in the inset is for an a-axis crystallite outgrowth, both showing a 
fairly low zero bias conductance. The data in[!J> was acquired on an a-outgrowth from a DC 
sputtered film. The relatively small gap in this figure may be due to local deviations from 
optimal doping resulting from losses of oxygen from the surface. The gap distribution in these 
regions were much larger as compared to the smooth ones, varying between 12 to 17 meV, 
with no apparent dependence on film preparation method. 

As mentioned earlier, the U-shaped gaps were found only on smooth areas, whereas V- 
shaped structures are typically measured on more disordered regions. A question then arises 
regarding the role played by the local surface morphology on the measured spectra. Since the 
typical gap size for both gap shapes (and also in the presence of a ZBCP, see below) is the same, 
around 17meV, it appears that the local surface roughness does not affect much the order 
parameter at the surface. We note in passing that our data are in agreement with Ref. |3K|. 
where similar gap values were obtained for both [110] and [100] orientations, further confirming 
the conclusion that YBCO is in the d-w&ve weak coupling regime. However, surface 
disorder is known to relax conservation rules for the in-plane component of the fc-vector, and 
in our case, may wash out the strong fc-dependent tunneling probabilities introduced in either 
the 'orbital coupling' or 'narrow tunneling cone' models. Surface roughness seems also to affect 
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Fig. 4 - dl/dV vs. V curves measured on a faceted region on the YBCO film shown in fig. The 
curves can be fitted well by a 'two channel' tunneling model, assuming independent contributions of 
tunneling in the [110] and [100] directions. The same A = 17 mV and V — 0.15A were used in all 4 
panels, and the different relative contribution of [110] tunneling is indicated in each panel. The fits 
presented here were calculated with the 'narrow tunneling cone' model, but similar fits are achieved 
with the 'orbital coupling' model. 



the quasi-particle life-time at the surface, smearing the gap structure and thus enhancing the 
Dynes parameters needed for the fits. 

Another typical spectral feature that we have found is a ZBCP embedded in a gap struc- 
ture, resembling dl/dV-V characteristics obtained for macroscopic tunnel junctions on (100) 
YBCO surfaces |13I24| . This behavior was attributed to the effect of (110) facets. Indeed, STM 
measurements have shown that even a facet of unit-cell height can give rise to a ZBCP |14I28| . 
Such a local sensitivity was not yet demonstrated for the (100) YBCO surface. In fig. 0]we 
present spectra exhibiting 'peak in a gap' structures measured on a faceted (rough) region 
of the laser ablated film shown in fig. ^i, where the relative ZBCP height increases from 
St (no visible peak) toQJl. The curves can be fit well by employing a 'two channel' model, 
where independent contributions of tunneling in the [110] and [100] directions are taken into 
account. The magnitude of the [110] tunneling component varies spatially, as is indicated in 
each figure, whereas the gap and the Dynes parameters were nearly the same, ~ 17meV and 
0.15 A, respectively. The fit to the U-shaped gap in fig. 0Ji requires a fc-selective tunneling 
process. However, due to the relatively large zero bias conductance level (as compared to 
other U-shaped gaps), a small [110] contribution had to be included; (the peak structure is 
washed out due the broadening parameter). 

In conclusion, our data provide clear evidence that the local tunneling spectra acquired for 
a-axis YBCO films do not always reflect the bulk superconductor DOS averaged isotropically 
over the Fermi surface. Tunneling characteristics exhibiting U-shaped gap structures, rather 
than the V-shapes, are often measured. The observation of U-shaped gaps does not depend 
on the sample preparation method, but is correlated with the degree of surface smoothness. 
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Two possible models that reproduce equally well these spectra were considered: a narrow 
tunneling cone around the [100] normal direction, and a strong overlap of the tip electronic 
states with the d-wave electronic orbitals of the Cu ions emerging out of the surface, both 
enhance tunneling preferentially along the anti-nodal directions. 
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